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Abstract 

Quantum chemical calculations including embedding, scalar relativistic, and dynamic elec- 
tron correlation effects on Cs2NaYCl6:(CeCl6)^~ embedded clusters predict (i) red shifts of the 
d{t2gY transition with pressure and (ii) bond length shortening upon / — > d{t2g) excita- 
tion. Both effects are found to be connected which suggests that new high pressure spectroscopic 
experiments could reveal the sign of the bond length change. 
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I. INTRODUCTION 



Recent ab initio theoretical studies of the structure and spectroscopy of octahedral com- 
plexes of lanthanide/actinide ions have shown that the bond length between the /-element 
and the surrounding ligands shortens upon /" f"'~^d{t2gY excitation. ^'^'^'"''^ This is 
so for different halide ligands (F, CI, Br), environments (crystals, liquid solutions, gas 
phase), oxidation states of the ion (III, IV), and number of / electrons (across the 
lanthanide/actinide series). ^'^'^'^'^'^ Yet, it contradicts the widespread assumption that 
jn _^ f^'^d} excitations lengthen the impurity-ligand bonds. ^'^'^'^'^ 

Usually, the materials involving /-element ion complexes are liquid or solid solutions and a 
direct experimental proof of the bond shortening could be obtained, in principle, by means of 
ground (/") and excited state {f"'^^d^) EXAFS measurements. In this respect, a theoretical 
study of (CeXg)^" (X= F, CI, Br) complexes in cubic elpasolites, liquid acetonitrile solution, 
and gas phase, has pointed out the chloride and bromide complexes in acetonitrile as good 
candidates for excited state EXAFS experiments because the negative bond length shifts are 
largest and the liquid medium could favor the experimental setting.^ However, excited state 
EXAFS measurements are extremely demanding because the difficulties inherent to EXAFS 
experiments and their interpretation are extended by the need of pumping the samples onto 
the excited f"'~^d^ levels for long periods of time and by the uncertainties on the actual 
excited state population attained. As a matter of fact, no successful experiments of this 
kind have yet been reported in octahedral complexes of / element ions in solid or liquid 
media, to our knowledge. 

In this paper, on the basis of quantum chemical simulations, we propose that spectroscopic 
studies under high hydrostatic pressure^^'^^ are an alternative, relatively simple experimental 
technique to proof or reject the predicted bond length shortening upon /" f"'~^d{t2gy 
excitation. In effect, quantum chemical calculations in Cs2NaYCl6:Ce^"'" reveal that the 
shortening in bond length upon /" f'^'^d{t2gY excitation results in a continuous redshift 
of the transition energy with pressure in the range 1 bar - 26 kbar. More generally, the 
sign of the bond length change upon excitation determines the sign of the slope of the 
transition energy with pressure. As a consequence, the theoretically found sequence of bond 
lengths in octahedral complexes.'^ R(\f^~^ d{t2gY] < Re[f"'] < Re[f"'~^d{egy], consistently 
leads to deacreasing /" f^'^ d{t2gY ■, constant /" — > /", and increasing /" f'"~^d{egy 
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transition energies with hydrostatic pressure. These results also enhance the value of high 
pressure experiments to alter the relative positions of excited states^'^^'^^ and to help the 
usually difficult assignment of the f^'^d^CgY highest energy levels. 

II. RESULTS AND DISCUSSION 

The effects of hydrostatic pressure on the local structure and electronic transitions of the 
(CeCle)^" defects in Cs2NaYCl6 were studied using the relativistic ah initio model potential 
(AIMP) embedded cluster method. ^^'^^ More detailed descriptions of the method, as apphed 
to / clement ions in ionic crystals, can be found in Refs. 2 and 4. The effective core po- 
tentials used for Ce and CI in an experimental and theoretical study of the absorption 
and emission spectra of Ce^"*" in elpasolite lattices^'' were also used here to incorporate scalar 
relativistic effects. In order to include dynamic electron correlation, SCF calculations were 
performed on the 4/^ and hd} electronic states and were used as reference for second order 
perturbation calculations where 57 valence electrons were correlated that occupy molecular 
orbitals with main character Ce 5s, 5p, 4//5(i and CI 3s, 3p.^^'^° (The SCF results of the 
structure and pressure-induced shifts of electronic transitions are in qualitative agreement 
with those presented here; therefore, the predictions made in this paper already hold quali- 
tatively at that level of theory.) The effects of pressure were modelled by using the AIMP 
embedding potentials produced in a recent study of high pressure effects on the structure and 
spectroscopy of defects in Cs2NaYCl6.^^ The embedding potentials corresponding to a 
given lattice constant were obtained through self-consistent embedded ions calculations;^^ 
they allow to incorporate quantum mechanical interactions with the host Cs2NaYCl6 ions 
in the (CeCle)'^^ cluster Hamiltonian.^^'^^ The values of the lattice constant used are listed 
in Table I, where the corresponding —l\V/V values are also shown. These data define a 
pressure range 1 bar - 26 kbar (see Table I) if the isothermal compressibihty of the material 
is K = (-AT//y)/AP = 4x10-3 kbar-^21 

Since Ce^"*" substitute for Y^"*" ions that occupy Oh sites the embedded cluster energies of 
all electronic states in the f"^ (the ground state '^A2u-, ^72«) and ^Ti^) and (^^2g and "^Eg) 
configurations were calculated at different Ce-Cl distances, R, for each lattice compression 
(there has been no reports of significant coupling to Jahn-Teller active normal modes). The 
plots of the potential energy surfaces for ambient pressure and 26 kbar are given in Fig. 1. 
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We obtained the equilibrium distances, Re, and totally symmetric vibrational frequencies, 
I'aig, as in Ref. 21; the vertical (AE') and minimum-to-minimum (AE'e) transition energies, 
and the shifts of bond length that occur upon excitation (Ai?e) were also obtained and are 
presented in Tabic 1 and Figures 1 to 3. 

The results in Table I and Fig. 1 show that the Ce— CI bond length shortens upon 
'^A2u — > d^ '^T2g excitation at ambient pressure and at higher pressures as well. Figures 1 
and 2 reveal a redshift of the ^^2^ — > d^ '^T2g transition with pressure. Both observations 
are connected: the sign of the bond length shift determines the sign of the slope of the 
electronic transition with pressure, being both of them negative, for the ^^2^ d^ '^'^2g 
transition. This connection can be established as follows. On the one hand (see Fig 3), the 
slope of the calculated transition energy AE with the Ce— CI distance R, dAE/dR, varies 
little with pressure (see solid lines) and has the same sign as (and is very similar to) the slope 
of the minimum-to-minimum transition energy AE^, with the Ce— CI equilibrium distance 
of the ground state, dAEg/ dRe(^ A2u) ■ The same is true for the slope of the vertical, Frank- 
Condon transition A£'[i?e(^^2«)] with Re(^A2u), dAE[Re(^A2u)]/dReCA2u)- Furthermore, 
the observations and conclusions drawn from Fig 3 also apply to any other i ^ f transition 
studied here, so that 

. dAEji ^ /) . dAEeji ^ /) 

dR = of , (1) 

and their magnitude is similar. On the other hand, it has been shown^ that when the force 
constant of the initial and final states are similar, ki ^ kf ^ k, like it is the case here (note 
that Paig values are given in Table I and k — Mci^-K'^Vaig^c^ for the aig normal mode) the 
slope of the transition energy dAE(i — > f)/dR depends basically on the force constant k 
and on the bond length shift ARg{i — > /) = Re{f) ~ Re{i)-i accompanying the transition, as 
follows:^ dAE{i f)/dR^ —kARe{i f) . This, together with the previous observations 
(Eq. 1), and the fact that dR^^i/dP < 0, leads to the conclusion that the sign of the slope of 
a transition with pressure coincides with the sign of the bond length shift that occurs upon 
excitation: 

sign of "^^^^^p^ = sign of Ai?,(i ^ /) , (2) 

and the same is true for the Frank-Condon transition. In effect (Fig 2, Table 1), the / — / 
transitions (with AR^ k, 0) are quite insensitive to pressure, / — > d{t2g) transitions (with 
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Ait!e < 0) shift to the red, and / — > d{eg) excitations (with AR^ > 0) shift to higher energies 
with pressure. 

It is known, from transition metal spectroscopy, that the so called lODg transition, 
d{t2g)^d{eg), shifts to high energy values with pressure, ^^'^"^ which is related to the bond 
length increase upon d{t2g)^d{eg) excitation (Eq. 2); the same behavior is obtained here 
for f''-'^d{t2gY f^-^diegY (Table I, Fig 2). It should be noticed, however, that this 
information alone does not suffice to infer the effects of pressure on the / — > d transitions; 
rather, the position of both potential energy surfaces f'^^^d{t2gY ^"^^ f"'~^d{egY relative to 
the potential energy surface of the /" ground state in the R axis of the configurational dia- 
gram, has to be also known. These relative positions have been found to be Re[f"'^^d{t2gY] 
< Re[f^] < Re[f^~^d{egY] by calculations here and elsewhere, as commented above,^ and 
have been interpreted as resulting from the following simple model of interactions: The inner 
lanthanide (Ln) (or actinide. An) /" open-shell electrons are shielded from the ligands by 
the outer Ln (An 6p^) closed-shell electrons, whose interactions with the hgands deter- 
mine the bond distance in states of /" configuration. Instead, upon f-^f^'^d^ excitation 
one electron has crossed the 5p^ {6p^) barrier and becomes exposed to covalent interactions 
with the ligands, at the same time that it leaves a 4/ (5/) hole behind, available for charge 
transfer from the ligands, both of which contribute to shortening the bond length between 
the f"~^d^ baricenter and the hgands. Finally, a large d{t2g)-d{eg) octahedral ligand spht- 
ting shortens the f'^~^d{t2gY bond lengths clearly below the /" ones (Ai?e <0) and raises 
the f"'^^d{egY ones clearly above them (Ai?e >0).^ These bond length shifts ultimately 
lead to the predicted red shift of f^^d{t2gY and increase of f^^d{egY transitions with 
pressure in Cs2NaYCl6:Ce^"'". Should the sequence of bond lengths be that traditionally 
assumed: -Re[/"] < Re[f^~^d{t2gY] < Re[f^~^d{egY], both transitions would increase with 
pressure. Whether the former or the latter sequence of bondlengths actually occurs could 
be revealed by performing high pressure spectroscopic experiments on this model system 
Cs2NaYCl6:Cc'^+ in the range 1 bar - 26 kbar. 

From all the results we have discussed here it is possible to derive a rule of thumb that 
can be applied to predict the sign of the change of a transition with pressure. It consists 
of using the calculated configurational diagram for ambient pressure (dashed lines in Fig 1, 
for instance) and use it assuming that the effects of pressure simply correspond to moving 
inwards across the R axis. This approximation applies as a consequence of Eq. (1), which 
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is a manifestation that the major effect of pressure on these local electronic transitions is a 
moderate decrease of the f-element - ligand distance. 

III. CONCLUSIONS 

Previous ab initio wavefunction based studies of the structure and spectroscopy of octahe- 
dral complexes of / element ions in crystals, liquid solutions, and gas phase have shown that 
the bond length between the /-element and the hgands shortens upon /" — > f'^~^d{t2gY 
excitation. This result contradicted the widespread assumption that the bond length is 
longer in the f^^^^d^ states. No experimental proof of the actual sign of the bond length 
change upon excitation has been given. We have presented here a quantum chemical simula- 
tion of spectroscopic experiments under high hydrostatic pressure. The simulation includes 
scalar relativistic effects, second order perturbation treatment of the valence electron corre- 
lation, and quantum mechanical embedding techniques. A continuous redshift of the lowest 
jn _^ f'^-'^d} band of Cs2NaYCl6:Ce^+ with pressure is predicted in the range 1 bar to 
26 kbar. The pressure-induced redshift is shown to be related to the bond length shortening 
that accompanies the / — >■ d{t2g) excitation. Therefore, the simulation points to the experi- 
mental detection of the pressure induced redshift as a way to support (or reject) the bond 
length shortening. More generally, the calculations suggest that the sign of the bond length 
change upon excitation determines the sign of the pressure-induced shift of the transition, 
and the effects of pressure on the / — > / (no shift) and / — > d{eg) (shift to higher energies) 
transitions are also predicted and are related to the respective negligible and positive changes 
in bond length. The Cs2NaYCl6:Ce^+ crystal is pointed out as an ideal model system for 
the counterpart high pressure experimental studies because the higher d{egY band 

has been observed below the cutoff of the host absorption. 
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TABLE I: Effects of high pressure on the potential energy surfaces of 4/^ and 5d} manifolds of 



Cs2NaYCl6:Ce^+. Distances in 


A, vibrational frequencies 


and energies in 


cm ^. 




lcittic6 coiistcint, (J 


10.7396 


10.6752 


10 5463 


10.4067 


10 3530 


-AV/V 





0.0179 


0.0530 


0.0902 


0.1042 


pressure 


1 bar 


4.5 kbar 


13.2 kbar 


22.5 kbar 


26 kbar 




2.687 


2.680 


2.664 


2.646 


2.638 




2.688 


2.681 


2.665 


2.647 


2.639 


2rp 

_/ 1 -1. 


2.690 


2.683 


2.667 


2.649 


2.641 




2.645 


2.639 


2.625 


2.609 


2.602 




2.705 


2.698 


2.682 


2.662 


2.655 




2.669 


2.663 


2.648 


2.630 


2.526 


Ai2p 4fi ^2,,^ 4fi 2T2,, 


0.001 


0.001 


0.001 


0.001 


0.001 


-L ^ti 

LU 


0.003 


0.003 


0.003 


0.003 


0.003 




-0.042 


-0.041 


-0.039 


-0.037 


-0.036 


y 


0.018 


0.018 


0.018 


0.016 


0.017 




-0.018 


-0.017 


-0.016 


-0.016 


-0.015 


22^2s^ ^Eg 


0.060 


0.059 


0.057 


0.053 


0.053 


dig J -^'■ZU 


306 


311 


322 


335 


341 


2t 

J. Oo. 


306 


311 


323 


336 


342 




307 


312 


323 


336 


342 




307 


312 


323 


336 


341 


'Eg 


300 


305 


317 


330 


336 


AEe 4fi ^^2,,^ 4fi 


460 


460 


480 


490 


500 


ILL 


890 


910 


950 


1000 


1020 




24300 


24100 


23700 


23200 


23000 


2F 


47200 


47300 


1 uu 




^oZUU 




33460 


33380 


33300 


33160 


33080 


22^2g^ Sd^ 2^;^ 


22900 


23200 


24000 


24900 


25200 



"Assuming {-/W/V)/P = 4 x lO'^ kbar"!. 

^Re of bd^ center of gravity = [3 x i?e(^?2g) + 2 x Rei^Eg)] /5 

<=A£;e of 4/1 ^ 5di center of gravity = [3 ^^E^{^A2u Tag) + 2 x AE^{^A2u Eg)] /5 



FIG. 1: Effects of pressure on the potential energy surfaces of 4/^ and 5d electronic states of 
Cs2NaYCl6:Ce^+. High pressure curves (solid lines) correspond to —AV/V = 0.1042 (26 kbar, if 
K = {-AV/V)/P = 4x10-3 kbar-i). 

FIG. 2: Effects of pressure on the 4/^— *^4/^ and Af^^5d} electronic transitions of 
Cs2NaYCl6:Ce3+. 

FIG. 3: Solid lines: 4/^ ^^2w— '^T2g transition energies versus Ce-Cl distance at different 
pressures (1 bar and 4.5, 13.2, 22.5, and 26 kbar, respectively, starting from the upper line). 
Frank-Condon transition energies, E[^T2g\Re(^ A2u)] — Ep742„;i?e(^^2«)] (circles) and minimum 
to minimum energies, E[^T2g;Re{^T2g)] — Ep^2u;-Re(^^2«)] (squares), versus i?e(^^2«) values at 
different pressures. 
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Figure 1. Ruiperez et al. Journal of Chemical Physics 
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Figure 2. Ruiperez, et al. Journal of Chemical Physics 
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Figure 3. Ruiperez, et al. Journal of Chemical Physics 
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